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(Dated: February 7, 2008) 

We search for the charmed pentaquark candidate reported by the HI collaboration, the Oc(3100) , 
in e^e~ interactions at a center-of-mass (cm.) energy of 10.58 GeV, using 124fb~'^of data recorded 
with the BABAR detector at the PEP-II e^e~ facility at SLAG. We find no evidence for such a state 
in the same pD*~ decay mode reported by HI, and we set limits on its production cross section 
times branching fraction into pD*~ as a function of cm. momentum. The corresponding limit on 
its total rate per e^e~ — > qq event, times branching fraction, is about three orders of magnitude 
lower than rates measured for the charmed Ac and Sc baryons in such events. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Ten experimental groups have recently reported nar- 
row enhancements near 1540 MeV/c^ in the invari- 
ant mass spectra for nK~^ or pK'^ jj]. The minimal 
quark content of a state that decays strongly to nA'+ 
is dduus; therefore, these mass peaks have been inter- 
preted as a possible pentaquark state, called 9(1540)+. 
The NA49 experiment has reported narrow enhance- 
ments near 1862 MeV/c^ in the invariant mass spec- 
tra for S'~7r~ and S'^7r+ |3; the former has minimal 
quark content dssdu, and these two mass peaks have also 
been interpreted as possible pentaquark states, named 
=■(1860)— and ^(1860)" [also known as $(1860)], with 
the latter being a mixture of ussuu and ussdd. The HI 
experiment has reported a narrow enhancement at a mass 
of 3099±6 MeV/c^ in the mass spectrum for pD*~ Q, 
which has a minimal quark content of uuddc, making this 
a possible charmed pentaquark state, named 8c(3100)". 
On the other hand, there are numerous experimental 
searches with negative results ]^: several experiments 
observe large samples of strange baryons with mass sim- 
ilar to that of the 9(1540)+, e.g. A{1520)^pK-, but no 
evidence for the 8(1540)+; several observe large samples 
of the nonexotic S~ baryon, but not the H'(1860)~~ or 
^(1860)° states; and several with large samples of D*~ 
do not observe the 8c(3100)° state. Our recent search j3] 
for the 6(1540)+ and S'(1860)~" in e+e~ annihilations 
found no evidence for these states, and we set limits on 
their production rates in e+e~ —^ qq events of factors 
of eight and four, respectively, below rates expected for 
ordinary baryons of the same masses. 

Here we report the results of an inclusive search for the 
charmed pentaquark candidate 8c(3100)'^ in e+e~ anni- 
hilation data; we expect equal production of the charge 
conjugate state, and its inclusion is implied throughout 
this article. The data were recorded with the BABAR de- 
tector ;.6j at the PEP-H asymmetric-energy e+e~ storage 
rings located at the Stanford Linear Accelerator Center. 
The data sample represents an integrated luminosity of 



124 fb collected at an 



cm. energy at or just be- 
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low the mass of the T(4S') resonance. We study the same 
decay mode as in the HI analysis, 8c(3100)° — > pD*~ , 

where the D*~ decays to D -kJ [ttJ denotes a "slow" 



W5 



pion from the Z?* decay), and the D decays to iiT+Tr . 

In addition, we consider the mode in which the D decays 
to i^+7r~7r+7r~. 

The BABAR detector is described in detail in Ref. |a. 
We use all events accepted by our trigger, which is more 
than 99% efRcient for both e+e~ -^ qq and e+e~ —> T(45') 
events. We use charged tracks reconstructed in the five- 
layer silicon vertex tracker (SVT) and the 40-layer drift 
chamber (DCH). The combined momentum resolution, 
ct(pt), is given by [a{pT)/pTf = [0.0013pt]^ + 0.0045^, 
where pt is the momentum transverse to the beam axis 
measured in GeV/c. Particles are identified as pions, 
kaons, or protons with a combination of the energy- 
loss measured in the two tracking detectors and the 
Cherenkov angles measured in the detector of internally 
reflected Cherenkov radiation (DIRC). 

We evaluate the 6c(3100)'' reconstruction efficiency 
and invariant mass resolution from two simulations. For 
production in e+e^— > cc events, we use the JETSET Q 
Monte Carlo generator with the mass and width of the 
i;c(2455)° baryon set to 3099 MeV/c^ and 1 MeV, re- 
spectively, and allow only the pD*^ decay mode. We 
leave all other parameters unchanged, and a momentum 
spectrum similar to those of nonexotic charmed baryons 
is produced. The events have a total charm of ±2, but 
this has negligible effect on the number and distribution 
of additional particles in the event, which are the quan- 
tities of interest here. We also simulate T{AS) decays in 
which one B decays generically in our standard frame- 
work [8J and the other decays into a state containing a 
Z'c(2455)'^ with parameters adjusted in the same way. 
This gives a much softer momentum spectrum, cut off 
at the kinematic limit for B meson decays, and a differ- 
ent environment in terms of other particles in the event. 
We find that the efficiency and resolution depend primar- 
ily on the 6c(3100)° momentum and polar angle in the 
laboratory frame, and negligibly on other aspects of the 
production process or event environment. We use large 
control samples of particles identified in the data to cor- 
rect small inaccuracies in the performance predicted by 
the GEANT-based Q detector simulation. 

We choose 8c(3100)'' candidate selection criteria de- 
signed for high efficiency and low bias against any pro- 
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FIG. 1: a) Distributions of the invariant mass difference be- 
tween the D*^ and D combinations in X^ pD*^ candi- 
dates in the data, for the two D decay modes, b) Distribu- 
tion of the reconstructed momentum of the D*~ in the e^e~ 
cm. frame for pD*~ candidates in the data (histogram); the 
dashed (dotted) line represents the D*~ spectra measured in 
T(45') decays (e+e~ -^ cc events) scaled as described in the 
text. 



duction mechanism. We use charged tracks reconstructed 
with at least twelve coordinates measured in the DCH, 
and select identified pions, kaons and protons. The iden- 
tification criteria for pions and kaons are fairly loose, 
having efhciencies better than 99% and misidentifica- 
tion rates below 1% for momenta below 0.5 GeV/c where 
energy loss in the SVT and DCH provide good separa- 
tion, and efficiencies of roughly 80% and misidentification 
rates below 10% for momenta above 0.8 GeV/c where the 
Cherenkov angles are measured well in the DIRC. The 
criteria for identified protons are tighter. For momenta 
below 1 GeV/c and above 1.5 GeV/c the efficiencies are 
better than 95% and 75%, and the misidentification rates 
are below 1% and 3%, respectively. 

In each event we consider every combination of identi- 
fied pii'~''7r~7r~ andpii'"''7r~7r^7r~7r~ and perform a topo- 
logical fit to each combination with the hypothesized de- 
cay chain X^pD*~'^pD 7r~ ^ pK~^Tr~ {tt^tt~)tt~ . No 
mass constraints are used in the fit, but the decay prod- 
ucts at each stage are required to originate at a single 
space point. The D has a finite flight distance, and we 
require the confidence level of the x^ for its decay vertex 
to exceed 10^^. 



We select candidates in which both the reconstructed 
D and D*~ masses are within 20 MeV/c^ of the peak 
value, namely 1843.8 < m;^+^-(^+^-) < 1883.8 MeV/c^ 
and 1989 < mj^+^_(^+^_)^- < 2029 MeV/c^. In Fig. ^ 
we show the distributions of the differences in recon- 
structed invariant mass Am = m^+^_^- — rrij^+T^- and 
m^f+Tr-TT+TT- for these X ^pK~^tt^tt^ 



'■K+ 



•K TT^TT TTs 



and pA'+TT^TT+TT^TrJ candidates, respectively. Clear sig- 
nals for D*^ are visible in both cases, with peak positions 
and widths (~0.6 MeV/c^) consistent with expectations 
from our simulation. The widths (~6 MeV/c^) of the 
corresponding D and D*~ peaks (not shown) are under- 
estimated by about 10% in the simulation. We require 
a mass difference within 2 MeV/c^ of the peak value. 



143.48<Am<147.48 MeV/c^. 

About 55,000 D*- -^ K+TT-TT- decays and 73,000 
D*~^ K~^7T~7T~TT^TT^ dccays are present in the selected 
data over respective backgrounds of 4,000 and 62,000 
random combinations. No event in either the data or 
simulation has more than one surviving pD*^ candidate. 
Without the proton requirement, over 750,000 D*~ are 
seen. Figure ^ shows the distribution of the D*~ mo- 
mentum, p*, in the cm. frame for the selected data. A 
characteristic two-peak structure is evident, in which the 
peak at lower p* values is due to D*^ from decays of B 
hadrons from T{AS) decays, and the peak at higher p* 
values is due to e+e~ — > cc events. For purposes of il- 
lustration, we show the spectra measured 10] from these 
two sources on Fig. ^, scaled by our integrated lumi- 
nosity, average efficiency and fraction of events with a 
proton. The shape is modified by the selection criteria; 
in particular, the proton requirement shifts the edge at 
the highest p* values. The background is verified by side- 
band studies to be concentrated at lower p* values; it is 
clear that we are sensitive to 8c(3100)° production from 
both of these sources. 

We evaluate the Oc(3100)° reconstruction efficiency for 
each search mode from the simulation, as a function of 
p* . High-mass particles at low p* are boosted forward 
in our laboratory frame, so that the probability of losing 
at least one track outside the acceptance is large, and 
the efficiencies are low, about 10% and 5% for the K^tt" 
and K^TT^TT^n^ modes, respectively. The efficiencies 
rise with increasing p* to respective maximum values of 
30% and 22% at the kinematic limit. The invariant mass 
requirements introduce negligible signal loss. The rela- 
tive systematic uncertainties on the tracking and particle 
identification efficiencies total 6-8%; at low and high p* 
values, there is a contribution of similar size from the 
statistics of the simulation. 

We calculate the 0c(31OO)° candidate invariant mass 
as mpD' = "^p/^+^-(^+^-)^J-"^K+^-(^+^-)^- + "^D-, 
where m]j*-= 2010 MeV/c^ is the known D*^ mass ^T]. 
We take the resolution on this quantity from the simula- 
tion, as it is insensitive to the simulated D'^*^ mass res- 
olution and previous studies involving protons combined 
with i^g [3 showed the proton contribution to be well 
simulated. We describe the resolution by a sum of two 
Gaussian functions with a common center. The width 
of the core (tail) Gaussian averages 2.5 (20) MeV/c^, al- 
most independent of p* , and the wider Gaussian con- 
tributes between 20% of the total at low p* and 10% at 
high p* . The overall resolution, defined as the FWHM 
of the resolution function divided by 2.355, averages 2.8 
and 3.0 MeV/c^ for the K~^7t~ and K'^7T~7r~^TT~ decay 
modes, respectively, with a small dependence on p* . 

We show nipD' distributions for the 0c(31OO)° candi- 
dates in the data in Fig. |2 for the two D decay modes. 
They show no narrow structure; in particular they are 
smooth in the region near 3100 MeV/c^, shown in the 
inset, where the bin size is two-thirds of the resolution. 
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FIG. 2: Invariant mass distributions for Oc(3100)'^ candi- 
dates in the data in the (black) K^tt^ and (gray) K^n~TT~-K^ 
decay modes, over a wide mass range and (inset) in the region 
near 3100 MeV/c^. 



Corresponding distributions for sidebands in the D and 
D*~ masses and the mass differences show overall struc- 
ture similar to that in the signal region. We consider 
several variations of the selection criteria that might en- 
hance a pentaquark signal, but in no case do we observe 
one. To enhance our sensitivity to any production mech- 
anism that gives a p* spectrum different from that of the 
background, we divide the data into nine p* ranges of 
width 500 MeV/c covering values from to 4.5 GeV/c. 
The background is lower at high p* , so we are more sen- 
sitive to mechanisms that produce harder spectra. There 
is no evidence of a pentaquark signal in any p* range. 

We quantify this null result by fitting a signal-plus- 
background function to the rUpD* distribution in each 
p* range. We use a p-wave Breit-Wigner lineshape 
convolved with the resolution function described above. 
The RMS width of the reported 9^(3100)° signal is 
12 MeWc^ and consistent with the HI detector reso- 
lution ^. Our mass resolution is considerably better, 
so we must consider a range of possible natural widths 
r of the 9c(3100)°. We quote results for two assumed 
widths, r = 1 MeV, corresponding to a very narrow 
state, and F = 28 MeV, corresponding to the width ob- 
served by HI, which we take as an upper limit. For the 
background we use the function ,f{m) — for m < ttiq 
and f{m) — -y/l — {mo/rn)^ exp(a[l — {mo / m)'^]) / m for 
m > niQ, where tuq — rup + mjj*- = 2948 MeV/c^ is the 
threshold value and a is a free parameter. We fit over 
the range from threshold to 3300 MeV/c^, except in the 
lowest p* range for the K^tt^tt^tt^ mode. Here the ac- 
ceptance drops sharply near threshold and the fit range 
is restricted to the region above 3000 MeV/c^. 

We perform maximum likelihood fits at several fixed 
9^(3100)° mass values in the range 3087-3111 MeV/c^. 
In every case we find good fit quality and a signal ampli- 
tude consistent with zero. We consider systematic effects 
in the fitting procedure by varying the signal and back- 
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ground functions and fit range; changes in the signal yield 
are negligible compared with the statistical uncertainties. 
The dependence on the assumed mass value is also small 
compared with the statistical error in each case. Fix- 
ing the mass to the reported value of 3099 MeV/c^, we 
obtain the event yields shown in Fig. |31 There is no pos- 
itive trend in the data, and the roughly symmetric scat- 
ter of the points about zero indicates little momentum- 
dependent bias in the background function. 

In each p* range we divide the sum of the two signal 
yields by the sum of the two products of reconstruction 



efficiency and D 



or D ^K^ 



branch- 



K+TT- 

ing fraction, the D*^ -^ D ■n'^ branching fraction, the 
integrated luminosity, and the p* range. This gives the 
product of the unknown 9c(3100)° -^ pD*~ branching 
fraction, B, and the differential production cross section, 
dcr/dp*. The resulting values oi B-Aa/dp* for r=l MeV 
and F=28 MeV are shown in Fig.0] We derive an upper 
limit on the value in each p* range under the assumption 
that it cannot be negative: a Gaussian function centered 
at the measured value with RMS equal to the total un- 
certainty is integrated from zero to infinity, and the point 
at which the integral reaches 95% of this total is taken as 
the limit. These 95% confidence level (CL) upper limits 
are also shown in Fig. ^ 

We integrate B ■ dcr/dp* over the full p* range from 
0-4.5 GeV/c, taking into account the correlation in the 
systematic uncertainty, to derive a total production cross 
section times branching fraction, B-a, for each of the two 
assumed F values, and calculate corresponding upper lim- 
its. These limits are model independent; any postulated 
production spectrum can be folded with the measured 
differential cross section to obtain a smaller limit. We cal- 
culate corresponding limits on the number of 90(3100)*^ 
produced per qq [q = udsc) event and per cc event by 
dividing by the respective cross sections for these types 
of events; we also calculate a limit per T(45') decay by 
integrating S- dcr/dp* over the range p* < 2 GeV/c (the 
kinematic limit for B meson decays is 1.8 GeV/c) and 
dividing by our effective cross section for e+e"— >T(4S'). 
These central values and limits are given in Table Q] 
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TABLE I: Total production cross section of the ec(3100)° 
pentaquark candidate times its branching fraction to pD*~ , 
Ba, in e'^e~ annihilations at ^/s =10.58 GeV, for two assumed 
values of the natural width. The corresponding 95% CL upper 
limits on B-a and on B times the yields per e"^e~ — > qq event, 
e^e~ -^ cc event, and T{4S) decay. 





r= 1 MeV 


r=28 MeV 


B-a (fb) 


40±44 
< 117 


102±111 
< 297 


S X yield X 10"'' per 
e e~ —> qq event 
e"'"e~— > cc event 
T{4S) decay 


< 3.4 

< 8.5 

< 12 


< 8.8 

< 22 

< 37 



In summary, we perform a search in e+e^ annihila- 
tions at -^/i =10.58 GeV for the pentaquark candidate 
state 9c(3100)° reported by the HI collaboration. We 



use the same decay mode as HI, 8c(3100)°^pl?*~, and 
find no evidence for the production of this state in a 
sample of over 125,000 pD*~ combinations. The com- 
ponents of this sample from c-quark fragmentation and 
B^ / B^ + B^ decays are both at least 100 times larger 
than the sample used by HI, implying that neither hard 
charm quarks nor B mesons produced in deep inelastic 
scattering can be the source of the HI signal. We set 
upper limits on the product of the inclusive 0c(31OO)'^ 
production cross section times branching fraction to this 
mode for two assumptions as to its natural width, which 
are valid for any state in the vicinity of 3100 MeV/c^. 
It would be interesting to compare these limits with the 
rate expected for an ordinary charmed baryon of mass 
^3100 MeV/c^. However rates have been measured for 
only two charmed baryons, the 71+ (2285) |10i, illj] and 
17c (2455) |l3i with precision that does not allow a mean- 
ingful estimate of the mass dependence. The mass depen- 
dence observed 11] for non-charmed baryons in e^e^ an- 
nihilations would predict a rate for a 3100 MeV/c^ baryon 
about 1,000 times smaller than that of the yl+(2285). 
Our limits for a narrow state in both e+e^ — > cc and 
T(4S') events are roughly 1,000 and 500 times below the 
measured yl+(2285) and i7c(2455) rates, respectively. As 
a result the existence of an ordinary charmed baryon with 
this mass and decay mode cannot be excluded. 
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